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Abstract 



■ - - Jladiative energy and momentum transfer due to fluctuations of electromagnetic fields arising due to temperature dif- 
ference between objects is described in terms of the cross-spectral densities of the electromagnetic fields. We derive 
relations between thermal non-equilibrium contributions to energy and momentum transfer and surface integrals of tan- 
gential components of the dyadic Green's functions of the vector Hclmholtz equation. The expressions derived here are 
applicable to objects of arbitrary shapes, dielectric functions, as well as magnetic permeabilities. For the case of radiative 
' transfer, we derive expressions for the generalized transmissivity and generalized conductance that are shown to obey 
["t I reciprocity and agree with theory of black body radiative transfer in the appropriate limit. 
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E Electric field vector 

tF View factor 

Ge Electric dyadic Green's function 

Grn Magnetic dyadic Green's function 

V X 

Gm V X Gm 

'Go Green's function of contribution due to back- 
ground or source radiation 

■=(sc) 

'G Green's function of contribution from waves scat- 
tered by interfaces 

G'' Linearized conductance for radiative transfer 

H Magnetic field vector 

/ Identity matrix 

J Current Density 

Ti Temperature in object I 

P Poynting vector 

Q Radiative heat transfer 

Rhn Fresnel reflection coefficients at interfaces between 
h and n 

Si Closed surface of object I 
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Generalized transmissivity for radiative energy 
transfer 

Generalized transmissivity for momentum transfer 
Volume of object I 

Volume of infinitesimal radius surrounding r 

Matrix of contribution to (EE*)^ 

Matrix of contribution to (HH*)s 

Matrix of contribution to (EH*)s 

Speed of light 

Reduced Planck's constant 

Wavevector 

Boltzmann's constant 

n component of wavevector (n = a;, y, z) 

z component of wavevector in vacuum 



1.2 



kl 



Thickness of vacuum gap 

Unit normal vector 

Position vector of observation point 

Position vector of source point 

Time 

or 1 

Delta function 
Distance of points 
Levi-Civita symbol 
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Tr 



Permittivity, 
Permittivity of free space 
Energy of a photon at temperature T 
Permeability, /i' + i/i" 
Permeability of free space 
1 or —1 

Maxwell stress tensor 
Frequency 
Real part 
Imaginary part 
Trace 



Superscripts 



hh 


Blackbody 


e 


Electric field 


m 


Magnetic field 


PP 


Planar- planar 


(h) 


Vacuum 


(0 


Objects (1 = 1,2,--- ,iV) 


ip) 


Transverse magnectic 


is) 


Transverse electric 


(m) 


Polarization s or p 


T 


Transpose 


* 


Complex conjugate 


Subscripts 



h Vacuum 

i,p, q Cartesian components 1, 2, 3 

/ Objects (? = 1,2,- •• ,iV) 

s Symmetric summation 

1 — > 2 From object 1 to 2 

1. Introduction 

Fluctuations of electromagnetic fields lead to thermal 
radiative transfer, via energy transfer, and van der Waals 
and Casimir forces, via momentum transfer. Diffraction 
and interference effects as well as tunneling of evanescent 
and surface waves, collectively known as near-field effects, 
are not taken into consideration by the classical theory 
of radiative transfer. Near-field effects become important 
when the length scale of importance becomes comparable 
to the characteristic thermal wavelength (A^ ~ SOOO/T 
fj,m). For radiative transfer between two objects, an impor- 
tant length scale is the minimum inter-object spacing, Igap- 
When Igap ^ At, tunneling of electromagnetic waves lead 
to enhancement of radiative transfer beyond the classical 
or far-field limit. Surface texturing, for instance by creat- 
ing a periodic ID or 2D pattern, introduces a length scale. 



Ip, that characterizes the period of the pattern. When 
Ip ^ At, diffraction effects can lead to thermal emission 
patterns not usually associated with a planar surface [ij. 

It has been long recognized that near-field enhancement 
of radiative transfer due to surface polaritons can result 
in increased power density as well as efficiency [l, H, 0, Hi • 
However, this enhancement of energy transfer has not been 
used in any practical device, as yet, because of our inability 
to conceive of configurations other than two parallel sur- 
faces with a thin vacuum gap in which an enhancement of 
similar magnitude occurs. Most investigations of near-field 
radiative transfer have been restricted to objects of few 
simple geometric shapes, each analyzed by a vector eigen- 
function expansion method applicable to that geom etry 
(planar geometry with vector plane waves 0, 0, mM, II, [lo| . 
cylindrical surfaces with vector cylindrical waves [ll| , two 
spheres with vector spherical waves [HI, [l^, H, , sphere- 
plane with a combination of vector spherical and plane 
waves [lij). Even minor changes to the shape of the ob- 
ject can impose great challenges. Simulations of thermal 
emission from textured surfaces are usually performed us- 
ing rigorous coupled wave analysis (RCWA) [T R [isl . Hifor 
finite difference time domain (FDTD) methods [20[, which 
are quite different from those used for simulations of near- 
field radiative transfer. To design other types of surfaces 
that can exploit the enhancement, without posing the hur- 
dles associated with two parallel surfaces, and also to de- 
sign surfaces with new radiative properties by shape mod- 
ification at nano/micro scale, we need a general method to 
predict all types of nanoscale effects on radiative transfer, 
irrespective of the size, shape or properties of the objects 
involved. 

Kruger ct al. [ll. 
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used fluctuational electrodynam- 
ics to develop a scattering matrix and operator formal- 
ism for computing non-equilibrium force and heat trans- 
fer interactions between objects with arbitrary shapes and 
frequency dependent dielectric permittivities. Biehs et. 
al. [9[ developed a formalism of nanoscale radiative trans- 
fer between two parallel surfaces similar to that of Lan- 
dauer formalism of electron transport in mcsoscopic de- 



vices 



25| 



26[ used 



Ben-Abdallah et. al. 
Rytov's theory to develop a theoretical formalism for ra- 
diative transfer between many objects in the dipole limit. 
Messina ct al. (2^ proposed a scattering matrix version 
of nanoscale radiative transfer as well as dispersion forces 
that is valid for objects with arbitrary shapes as well as 
dielectric functions. Non-equilibrium fluctuational electro- 
dynamical interactions between objects can be expressed 
in a scattering matrix formalism or in a Green's functions 
formalism, just as the electrical conductance for electron 
transport can be developed in terms of the scattering ma- 
trix or Green's function. 

The work in this paper is an extension to a prior work 
published in this journal by one of the authors [28| . In Ref. 
[28|, the focus was on the relation between cross-spectral 
densities of electromagnetic fields in thermal equilibrium 
and the dyadic Green's functions (DGFs) of the vector 
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Figure 1: Schematic of N objects at temperatures Ti,T2,--- ,Tjq 
embedded in a host medium at T}^. 



Helmholtz equation. In this paper, the focus is on ther- 
mal non-equilibrium effects, i.e. when the objects are at 
different temperatures. Volume integral expressions for 
cross-spectral densities of components of the electric and 
magnetic fields that can be obtained from Rytov's the- 
ory of fluctuational electrodynamics are converted into a 
form more appropriate (in terms of surface integrals of 
DGFs of the vector Helmholtz equation) for computations 
as well as comparison with the classical theory of radia- 
tive transfer. Though the focus of this paper is not on 
developing new numerical techniques, it is hoped that the 
formalism developed here will be used to compute thermal 
non-equilibrium energy and momentum transfer between 
arbitrarily shaped objects. 

The paper is arranged as follows. In Sec. [51 the 
fluctuation-dissipation theorem and DGFs are used to ex- 
press the electric and magnetic field correlation functions 
in terms of the volume integrals of the DGFs. In Sec. [31 
Green's identities for dyadic functions arc used to derive 
expressions for the field correlations in terms of surface in- 
tegrals of tangential components of the DGFs. For radia- 
tive transfer between two objects, a generalized transmis- 
sivity function that is expressed in terms of double surface 
integrals on the surfaces of the two objects is derived in 
Sec. [H In Sec. [3 the theoretical formalism developed in 
Sec. [H is applied to the cases of radiative heat transfer 
between two parallel half spaces. We also show that the 
generalized transmissivity function agrees with the theory 
of blackbody radiation in the appropriate limit. Finally, 
we also discuss the implications of the theoretical formal- 
ism developed here for computation of heat transfer and 
non-equilibrium forces. 

2. Fluctuational electrodynamics and Green's 
function formalism 

We briefly describe our notation regarding electromag- 
netic fields and their Fourier transforms here. A field 
Aiv^t) and its Fourier transform, A(r, cj), are related by 



is used to identify a field as well as its Fourier transform, 
explicit dependence on time will be included when rcfering 
to the time domain field. Explicit dependence on lu is sup- 
pressed from Aiv^Lo) so that it is written as ^(r). Explicit 
dependence of relative dielectric permittivities, magnetic 
permeabilities, and DGFs on is also suppressed. 

Let us consider N objects (see Fig. [Ij with relative 
dielectric permittivities si{uj) and magnetic permeabilities 
/i;(aj) at temperatures T;, where I = 1,2,-- - ,N. These 
objects are assumed to be embedded in vacuum that is at 
temperature T/j. The object I is confined to the volume 
Vi and the closed surface Si is the boundary of this object 
with the host medium. The outward normal on the surface 
of object / at r is represented by n; (r) . The fluctuations of 
the electric and magnetic current densities, which give rise 
to the dispersion forces and radiative transfer, are related 
to temperature by fluctuation-dissipation theorems of the 
second kind (29l . 30l . 3l|: 



(j;"(r)J™*(r)) 
(j;(r)J™*(r)) 



2a;eo£"e(w,r)(5(r-r)(5pg (la) 

2LJ^lo^Ji"e (w, t) 5{y - v)5pq (ib) 

(Ic) 



where p^q = 1,2,3 are the labels for the Cartesian com- 
ponents of the vector, Eq and /io are the permittivity and 
permeability of free space, and J™ are the Cartesian 
components of the electric and magnetic current densities, 

hw f huj \ 
S (uj,T) = — coth , e and /.i are the imagmary 

2 \2kbi J 

parts of the dielectric permittivity and magnetic perme- 
ability respectively at the location r which is in local ther- 
modynamic equilibrium at temperature T, z* is the com- 
plex conjugate of z, and () denotes the ensemble average. 
27rft is the Planck constant and fcf, is the Boltzmann con- 
stant. The presence of Spq implies that we assume all ma- 
terials to be isotropic, and that of 5(r — r) implies that the 
correlations of sources are local. The Fourier transforms of 
the electric and magnetic fields in the host medium (vol- 
ume Vh in Fig. [IJ due to sources in object I (volume Vi) 
are given by: 

E(r) = J [p(r) . ^e(r, r) - J™(r) • ^^(r, r)]dr (2a) 

Vi 

H(r) = I [m(r) • ^„,(r, r) + J^(r) • ^A,(r, ? )]dr (2b) 

Vi 

where p {r)~iuj^ofJ-ir)J^ (r), m (r) = iujeo£{r)J"^ (r), r G 
Vh,GEir,r) ^VxGe(r,f) andGA/(r,f) = VxG™(r,f). 
Ge(r,f) and Gm(r,f) are DGFs of the vector Helmholtz 
equation that satisfy the following boundary conditions on 
the interface Si between object I and the host medium: 

n;(r;) X (^;(r/)Ge(r/,r) - ^ft(rft)Ge(r,i, r)) = 0, (3a) 
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hii^i) X (Gi;(rz,f)-G£;(r;„r)) = 0, (3b) 



hi(ri) X (e/(ri)G,„(r/,r) - e,,(r/,)G„(r,i,r)) = 0, (3c) 



M^i) X (GM(r,,r) - GM(r,„r)) = 0, (3d) 

where r; and Vh are position vectors of points on either 
side of 5; in volume Vi and Vh respectively (jr; — r^| — >■ 
0) . In addition, the DGFs satisfy the following reciprocity 
relations: 



n{r)G^ (r, r) = /i(r)Ge(r, r), 



r(r)G™(f,r), 



G'£;(r,f) = G'M(f,r), 



(4) 



(5) 



(6) 



where A is the transpose of A. 

Radiative transfer can be determined from the Poynting 
vector, P(f) = (E(r,t) x H(r,t)), whose components are 
given by: 

P,ir) ^ e,p,{Epir,t)H,{r,t)) (7) 

where eipq is the Levi-Civita symbol. To determine 
van der Waals pressure and radiative transfer, we need 
equal time correlations of various components of the elec- 
tric and magnetic field vectors, such as {Ep{r,t)Eq{r,t)), 
{Hp{r,t)Hq{r,t)), and {Ep{r,t)Hq{r,t)). van der Waals 
pressure in vacuum can be determined from the Maxwell 
stress tensor, a = a +a , where ct and ct are the elec- 
tric and magnetic field contributions respectively, a and 
a are given by: 



r(f) 



cr (r) = /io 



(E(r,t)E(r,i))--/(E2(r,i)) 



(H(r,i)H(r,<))-i/(H2(r,t)) 



(8) 



(9) 



where (E(r, t)E(r, t)) and (H(f , i)H(f , t)) are ma- 
trices whose components are {Ep{r,t)Eq{r,t)) and 

{Hp{r,t)Hq{r,t)) respectively, p,q = 1,2,3, and / is the 
identity matrix. 

Since the fields are assumed to be stationary, 
{E(r,t)E(f,t)), (H(rJ)H(r,i)), and (E(r, t)H(r, t)) are 
independent of time |32| . The equal time correlation func- 
tions are related to the cross-spectral densities by: 



{Ep{r,t)Eq{v,t)) 



dio 
2^ 



{Ep{v)El{v)), 



(10) 



duj 
2^ 



{Ep{v)E;{v) + E;{v)Eq{v)) 



oo 

{Hp{v,t)Hq{v,t))^ J '^{Hp{v)H;{v)), 



(11) 



2^ 



Hp{v)H;{v) + H;{v)Hq{v)) 



{Ep{v,t)Hq{v,t)) 



2^ 



{Ep{i)H;{v)), 



OO 

:j ^{Ep{r)H;iv) + E;iv)Hqir)) 



(12) 



where the subscript s in Eq. [TU] - Eq. [T^ implies a sym- 
metric sum. Using Eq. [T] and Eq. ^ we can express the 
cross-spectral densities of the components of the electric 
and magnetic field at r e Vh as: 



(E(r)E*(r)), = 2c^A*o 
(H(r)H*(r)), = 2t^£„ 



N 



. , ={l) =(h) 

Y,Qi£"{v) + QhS (?) 



L;=i 

" N 



, , =(;) =(h) 

Y^em (f) + e„H (?) 



.1=1 



(13) 
(14) 



N 



(E(?)H*(r)) = J2 ^i^^'^ (f) + ©'^-^^'^ (f) (15) 

1=1 

where 6/ = e{u;,Ti), and (E(?)E*(r))„ (H(r)H*(r))„ 
and (E(r)H*(?)) are matrices whose components are 
{Ep{r)E*{r))s, {Hp{r)H*{r))s, and {Ep{r)H*{r)) respec- 
tivcly. Even though components of (EH*)^ arc necessary 
to compute radiative transfer, we persist with (EH*). The 
reason for computing (EH*) as opposed to (EH*)., will be 

clarified in Sec. 121 

=(0 =(0 =(i) 
The matrices £ (r), H (?), and X (?) are contri- 
butions to (E(?)E*(?)),, (H(?)H*(?)),, and (E(?)H*(?)) 
{I = 1, 2, • • ■ , A^, ft,) from sources in volume V/. For I G 

=(0 =(0 =(0 
{1, 2, • ■ • , N}, £ (?), n (?), and X (?) are given by: 

t\r) = 25R / dr\fi"{r0l{r,r)-%{r,r)+ 



(16a) 



e"irMr)\' — G, (r, ?) • GJr, ?) 

=(0 r r =T =* 

n (?) = 23? / dr e"(r)GM(r,?) • Gjv,(r,?)+ 



(16b) 



M"(r)|e(r)|2i^Gl(r,?).G„(r,?) 



x'^'^ (?) = ^2'^ I £"(r)/i(r)G, (r,?) • GM(r,?)+ 



(16c) 



M"(r)e*(r)Gs(r,?).G„(r,?) 



dr 



=(/i) =ih) =(h) 

The expressions for £ (r), (r), and X (?) have 
to be modified to take into account the singularity of the 
DGFs in the integrals in Eq. [TTallTTcl as |r - r| 0. The 

=(h) =(h) ={h) 

modified expressions for £ {v), % (r), and X (f) 
are: 



={h) 

£ (f ) = lim 25R 



M"(r)Gs(r,r).G^(r,f) 



-e"(r)|M(r)p^^C(r,f)-^I(r,f) 



dv 



(17a) 



+e 



=(sc)T , 



L-G, (?,?) + ^G, (r,r).L 



={h) 
n (r) 



lim 2» 

Vs^O 



e"(r)Gl,(r,r).GL(r,f) 



/i"(r)|£(r)|^ — G„(r,r).G„(r,r 



dr 



£f, =isc)T =- 

%G,„ (r,r).L 



(17b) 



=(/i) 
X (f ) 



i2 



a; 



lim 



e"(r)Mr)G, (r, r) • G^,(r, r) 



//"(r)e*(r)Gs(r,r).G,„(r,r) 



dr 



-z2(^L.G,, 



(f,f) • L 



(17c) 



where V5 is volume of infinitesimal radius surrounding 



r, L is a shape dependent dyad j33|, and Ge(r, r) = 
Go(r,r) + cl'^'^r,?), G„(r,r) = Go(r, r) + cj^'^r, r), 
Gi5(r,f) = V X Go(r,r) + G^ (r,r), and GM(r,r) = 
V X Go(r,f) + G]^j (r,r). Go(r,r) is the DGF when no 

={sc) 

scatterers are present, and G (r,f) is the contribution 

=(sc) 

from presence of scatterers. G (r, r) is always finite. 
The volume integrals in Eq. I17al -Eq. I17cl are finite even 
though Eft = = because of the singularity in the 
DGFs. The last line of Eq. HMEq. [T7c] are identically 
equal to zero for non-absorbing materials, including vac- 
uum. However, we choose to retain them since they are 
essential for calculations of cross-spectral densities in ab- 
sorbing media. We have neglected terms independent of 
the configuration of the scatterers. These terms are infi- 
nite because of the assumption that the thermal sources at 
any two locations are uncorrelated. They can be made fi- 
nite by eliminating the local assumption in Eq. [T]but that 
will not affect the calculations of forces or heat transfer ex- 
cept at gaps smaller than the correlation length. Usually, 
the correlation length is of the order of the atomic spac- 
ing in dielectrics or the electron mean free path in metals. 
More detailed discussion of the singularity in DGFs and 
calculation of cross-spectral densities is given in Ref. . 



3. Surface integral dyadic Green's function formal- 
ism 

While the volume integrals in Eq. II Gal - Eq. I17cl can 
in principle be used to compute forces and radiative heat 
transfer, they are undesirable for the following reasons: (1) 
Evaluating classical radiative transfer between two objects 
requires the computation of the view factor between them. 
But for objects with simple geometries, computation of the 
view factor between two objects requires, in general, not 
the evaluation of a volume integral but the evaluation of a 
double integral over the surfaces of the two objects, (2) the 
expressions in Eq. I16al - Eq. ll7cl do not reflect the different 
reciprocity relations and boundary conditions satisfied by 
the DGFs. and (3) evaluation of volume integrals are com- 
putationally more expensive than that of surface integrals. 
These undesirable features can be overcome by converting 
Eq. I16al - Eq. ll7cl into appropriate surface integrals using 
Green's theorems for dyadic functions [sj]- The surface 
integral representations for the cross-spectral densities are 
as follows: 

=(l) C =T = 

£ (r) = 23 i [^l{v)G^ (r,f)]*.[n,(r)xGi5(r,r)]dr (18a) 

Si 

=(i) r T 

n (?) =23 y [e(r)G„(r,r)]* • [n,(r) x GA/(r,r)]dr 
Si 

(18b) 



=(0 
X (r) 



Ge (r,r).(n,(r) xGM(r,r))- 



(18c) 



={h) 

£ (f) 



^^/i(r)Gr(r,f).(nz(r) x £(r)G„(r, ?))* 

N 



dr 



1=1 



2^3 f M*(r)Ge (r,r) • (nj(r) x G£(r,r))rfr 



+^i2nhGe{r,r)) 



(18d) 



={h) 
H (?) 



N 



1=1 



2^3 f £*(r)G,„ (r,r) ■ (nz(r) x GM(r,f))dr 



Si 



-f5(2e,,G,„(r,r)) 



(18e) 



N 

^^'''(r) pEir, r) ■ (n,(r) x %,{r, r))- 



'=^Si 



— M(r)G, (r,f).(nz(r) xe(r)G™(r,r))* 
i23(^M(f,f)) 



dr 



(18f) 



If the host medium is dissipative, the functions 
^{2^lh^,{r,r)), ^{2ehd„,{r,r)), and ^(^^(f,?)) in Eq. 
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Il8d[ Eq. I18el and Eq. Il8fl respectively should be replaced 

by 3(2/x,^l"\r,r)), 5(2£„^l7\f , ? )), ^nd^0lj \r,r)). 
Using Eq. llSalllSfl Eq. [HMTCl ean be re- written as: 



N 



object indices 1 and 2 can be replaced by any 
m,n G {1,2,- ■■ ,N}. From Sec. ^ we know that 

CJO 

dio 
2^' 



=(0, 



(£;,(?, t)ffp(r,i))^^^ = / — (i?,(r,c.)i7;(r,a.))W. Using 

1 Eq. m for Qi^2 



(EE*), =2wMo J2 " + 2w/ioe,,9(2/i,,Ge(r, r)) Eq. [13 Eq. HH Eq. W and Eq, 

z=i can be re-written as: 

(19) 

oo 

did 



N 



(HH*), 



=(0, 



(20) 



Qi^2 = - 



2^ e(w, Ti) ^ (irn2.(r)e,p,25RA'pW(f) 

S2 



( = 1 



(EH*) 



^(6, - eh)X^'\r) - i2e„3GM(f , r) (21) 



1=1 



The terms 5(2/i/iGe(f , f )) and SJ(2e;iGm(f , r)) are ther- 

={h) _ =(/») _ 
mal equilibrium contributions to £ (f) and Ti (f ) re- 
spectively. They give rise to the van der Waals stresses as 
predicted by Lifshitz theory when Sh = iJ^h = ^- The rea- 
son for persisting with (EH*), as opposed to (EH*)s = 
EH* + E*H, is to show that there is indeed an equilib- 
rium contribution to (EH*). However, when we compute 
(EH*)s, the equilibrium contribution vanishes since ra- 
diative energy transfer between two objects at the same 
temperature must be zero. 

Radiative transfer between two objects is discussed fur- 
ther in Sec. |H Before proceeding to Sec. HI we wish 
to remark on the form of the non-equilibrium contribu- 
tions in Eq. llSalllStl Using the property that A = 
n(^n- — n X n X ^, we see that all the surface inte- 
grals in Eq. llSalllStl feature only tangential components 
of the dyadic Green's functions that are continuous across 
an interface between two materials. 



4. Generalized transmissivity for radiative energy 
transfer 

The steady state radiative heat transfer from object 1 
to object 2 in Fig. [TJ Qi^2, is given by: 



Qi 



^P(^'(f) • n2(r)dr (22) 



where P'^^^(f) is the Foynting vector at r G 5*2 due to 
thermally fluctuating sources within Vi. The "— " sign in 
front of the surface integral is because n2 (?) is the outward 
pointing normal on the surface 82- The net heat transfer 
between objects 1 and 2 is given by Qi^2 = Qi^2 — Q2^i- 
The components of P(^^(r) are given by: 



,(i?,(r,t)iJp(r,t))(i), 



(23) 



where (i?q(r, t)ifp(r, t))'^^) is the contribution to 



{Eq{r,t)Hp{r,t)) from sources within Vi. 



The 



00 





>2(w) ^ - j dvn2^{v)e,r,q2V.Xll\T) 

S2 



-^dre,,pn2.(f)25RA'^i'(f) 
S2 

n2(r) X X (r) 



(24) 



-23? j,dv 

Si 



-25Rrr 



=(1) 

dv ( n2(r) X X (f^ 



S2 



where 7\%2('^) is a generalized transmissivity for radiative 

energy transport between objects 1 and 2, and Tr{A) = 
3 

^ App. The superscript e in Tf_j.2(a;) stands for "energy." 
p=i 

=(1) TTTT-i- 

Substituting the expression for X from Eq. I18cl in the 
last line of Eq. [211 Tf_^2(^) can be shown to be: 



(j) dr (j) 

Si S2 



dr 



U! 



[n2(r) X ^2Ge(f,r)] ■ [ni(r) x £*G,„(r,r)] 
+ [fi2(r) x G£;(r,r)] • [ni(r) x G^(r,r)] 



(25) 



where ni(r) is the outward pointing normal on the surface 
5*1 , as shown in Fig. [TJ For any two vectors a, b and 
dyads A, B, the following property can be shown to be 

true: Tr{ (a x A) • (b x B) } = Tr{ (b x A ) • (a x B ) }. 
Using this property, and the reciprocity relations (Eq. [31 
Eq. \5\ and Eq. [5]), we can derive the following equations: 



3fiTr[ri2(r) x M2Ge(r,r)] • [ni(r) x eJG„(r,r)] = 
JRTr [fii(r) x ^l{^e{r,r)] ■ [M^) x e2^lir,r)] 

3fiTr[n2(r) x ^i;(r,r)] • [ni(r) x ^^r, r)] = 
5Rrr[ni(r) x ^A/(r,r)] • [n2(r) x %j{r,r)] 



(26a) 



(26b) 



With the aid of Eq. [26l Tf_^2('^) can also be shown to be: 



" (j) dr (j) dr 

Si S2 



[rii(r) X ^iGe(r,r)] • [n2(r) x e^G„(r,r)] 
+ [fii(r) X GM(r,r)] • [n2(r) x GM(f,r)] 



(27) 



The generahzed transmissivity from object 2 to objeet 1, 
T2^-^{u!) can be determined from Eq. [57] (or Eq. ESI) by 
interchanging the subscripts 1 and 2 (r G S'l and r G S2 
are dummy variables and do not affect the value of the 
double integral). T|_^j(cj) is given by: 

^{uj) = 2'SiTr j dr j dv 



Si 



[fi2(r) X ^2Ge(r,r)] • [ni(r) x eJG„(r,r)] 
+ [n2(f) x GM(f,r)] • [ni(r) x G^^(r,f)] 



(28) 



That the expressions for generalized transmissivity de- 
rived earlier (Eq. [55] or Eq. 157)) satisfy the principle of 
reciprocity in thermal radiative transfer, i.e. r|_j]^(a;) = 
TiLj.2(a;), can be established by using Eq. I26bl to modify 
the expression for r|_^j(aj) as follows: 



T2^ii^) = 25ftrr (j) dr (j) dr 



S2 Si 

[112 (f) x /i2Ge(f,r)] • [ni(r) x £*G,„(r,r)] (29) 
+ [ni(r)x^£(r,r)].[n2(r)x^;(r,r)] 
= ri%2(^) (Eq. [23 

Though expressions for Tf_j.2(aj) in Eq. [55] and Eq. [57] 
are surface integrals, they are in fact derived from a volu- 
metric integral over Vi (Eq. I16cp . Similarly, energy emis- 
sion from the object Vi {I = 1,2, ■ ■ ■ , N) is derived from 
a volumetric integration over Vi. For this reason, the for- 
mulae derived for Ti%2('^) and Ti^i{uj) (Eq. [25l Eq. [ST] 
or Eq. [5^ can be described as "interior formulae." The 
integration over Vi (for T'f_^2('*^)) oi" ^2 (for T2_^i{uj)) is 
made explicit by the presence of £1, /i2 (in Eq. [55]) or £2, /^i 
(Eq. [57] Eq. [5^]). The corresponding "exterior formula" 
should not involve, or appear not to involve, any of these 
properties in the formula for transmissivity. The exterior 
formula for T^%2 can be derived by using the boundary 



conditions (Eq. [5alEq. 
the following equation: 



and converting Eq. [57] into 



rf_^2(^) = 2KT^ 



r <j) dr (j) dr 

Si S2 



[ni(r) X /iftGe(r,r)] • [112 (?) x £,;G,„(r,r)] 
+ [fii(r) X GA/(r,r)] • [n2(r) x Gj,f(f,r 



(30) 



There is a correspondence between the "direct" and "indi- 
rect" methods [35], [36], |37| and the "interior formula" and 
"exterior formula" derived above. It can be shown that 
the "exterior formula" is a generalization of the "indirect" 
method to include problems of near-field thermal radiative 
energy transfer between two objects, in addition to the cal- 
culation of thermal emission from objects for which it is 
currently used. This correspondence will be undertaken in 
a future work and is not pursued any further in this paper. 

Since r|_^]^(u;) = Tf_^2('^): the net radiative exchange 
between objects 1 and 2, Qi^2, is given by: 



n,2- / ^[e(a.,ri)-e(c.,T2)]Ti%2(^) (3i) 



From Eq. [JU a linearized conductance for radiative trans- 
fer between objects 1 and 2 can be defined as: 

00 

GUiT)= lim -91:^ = f ^^T,%2i^) (32) 



Generalized transmissivity for radiative momentum 
transfer? 

Though expressions for generalized transmissivity in 
terms of DGFs (Eq. [55] Eq. [57] Eq. [53 and Eq. ^ 
have been derived for energy transfer, we have been unable 
to obtain equivalent expressions for generalized (vectorial) 
transmissivity or conductance of thermal non-equilibrium 
momentum transfer. Why this is so can be explained by 
considering the nature of the Poynting vector and the elec- 
tromagnetic stress tensor. It is a well-known property 
of Maxwell's equations that the electric field at any lo- 
cation due to sources within a particular object, for in- 
stance object 1 in Fig. [TJ can be expressed in terms of 
surface integrals of tangential electric and magnetic fields 
on the surface of that object. This property of electro- 
magnetic fields forms the basis for the boundary element 
method for numerical solution of electromagnetic scatter- 
ing problems. The normal component of the Poynting vec- 
tor on the surface of object 2, (E(r, t) x H(r,i)) • n2(r), 
has an additional property that it can be written as 
[(n2(r) x E{r,t)) x (n2(r) x H(f , i))] • n2(r). This ensures 
that the radiative heat transfer between objects 1 and 2 
can be expressed in terms of tangential electric and mag- 
netic fields on the surfaces of both objects. This property 
of radiative heat transfer is reflected in the different ex- 
pressions for rf_j.2(w) because they contain only tangential 
components of the DGFs on the surfaces of both objects. 
However, the electric and magnetic stress tensors in Eq. 
[5] and Eq. [H] do not share this property. The force ex- 
erted by object 1 on object 2 requires knowledge of not 
only the tangential components of E and H on the surface 
of object 2 but also the normal components. While our 
inability to deduce an appropriate form for the transmis- 
sivity for momentum transfer does not mean that such a 
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transmissivity does not exist, Eq. [191 and Eq. [20] can still 
be used to determine ^'^(r) (Eq. [5]) and ^"'(r) (Eq. ^ 
for specific geometric configuration of objects, from which 
thermal non-equilibrium van der Waals forces between ob- 
jects can be computed. 



and 



riii-'-) — 



^hl ^h2 ^ 



if z^ = -l,C = l 



(36) 



5. Application to specific geometries or properties The total DGF, Ge (r, f), is given by Go (r, r) 4- G^, \r,r). 

The reflection coefficients and at the interface 



5.1. Planar multilayered media 

Since theoretical analysis of radiative transfer and ther- 
mal non-equilibrium van der Waals forces between planar 
multilayered objects have been published in literature, the 
expressions derived in Sec. [3] are applied to objects shown 
in Fig. [5] and the resultant expressions compared with 
those in literature. The half spaces in Fig. [21 marked 
1 and 2, can be homogeneous materials or can be com- 
posed of planar multilayer films. The only requirement is 
that the temperature gradients within 1 and 2 are negli- 
gible enough that they can be approximated as thermal 
reservoirs at temperatures Ti and T2. The vacuum layer 
of thickness / separating the two objects is equivalent to 
the host medium in Fig. [H The planar media lie in the 
X — y plane and the unit vector in the z direction, z, is 
directed from object 1 to object 2. The interfaces of ob- 
jects 1 and 2 with vacuum are at coordinates z = zi and 
z = Z2 {\z2 — zi\ = I). The polarization dependent reflec- 
tion coefficient of electromagnetic plane wave originating 
in vacuum and incident at the surface of half space 1 (in 
the absence of half space 2) is denoted as where 
fj, = s (transverse electric), p (transverse magnetic). A 
similar reflection coefficient for waves incident on half 
space 2, in the absence of half space 1, is denoted R'fl^2 ■ 
Go (r,r) for any two locations r and f within the vacuum 
layer is given by: 



Go(r,r) 



47r J khz 

x^>'H+khz)x^^'H+khz)e'^''~'''^'"'^ if z > z' 
i (-fc/iz)i*^' (-fc,iz)e^(^'-^)'='^- if z < z' 

(33) 



where kp = kpkp = k^x + kyfj, x^'^\±khz) ~ kp x z = 
{kyX-kxy)/kp, x'^P\±khz) = {Tkhzkp+kpz)/kh, and kl^ + 
kp = kf^. For vacuum, /i/i = e/j = 1 and kh — uj/c. The 

=(sc) 

scattered DGF, Gg (r, r), is given by: 



47r J khz 



e f 



G; '(r,r) - — / y -^x 



i/=±i4=±i 



where 



(34) 
(35) 



between two homogeneous media are given by the usual 
Fresnel reflection coefficients. For multilayered media, 
they can be computed using the transfer matrix method 
or by using recursion relations [s^. Since both r and r 
are within the same layer, G„i(r, r) and Gj\/(r,r) can be 
determined by simply replacing all occurrences of e(w) by 
the corresponding /u(a;) and vice versa. Explicit expres- 

sions for G^ (r,?), (r,?), and Gj^j (r, ?) are given 
below: 



=(sc) 

G^ (r,r) 



dk p k p 



khz 



E 



oikp-(p-p') 



E E Cl';^''S^''''H'^khz)x'^^\Ckh.)e'^-'-^''^''% 

(37) 

where fi' ~ p ii fi ~ s and fi' ~ s ii fi ~ p. Defining 
{kpkp ± khzZ)/kh, we have the following relations: 



k 



(±) 



(±) X ^(s) 



i±khz) 



hz , 



and Vf^ X (±fc. 



hz] 



{±khz). Using these relations, we obtain: 



Ge (r,f) = 



-1 
An 



dkpkp 
' kh 
khz 



E 



DM 



^ ^ (-l)^G(^^)£(^')(i^fc,,)5;(^)(efc,..)e'(''^-«^'"=- 

u=±l ^=±1 

(38) 



=(sc) 

Gm (r,i=) 



-1 f dkpkp ^ gik,-ip~p') 
4^7 "h^ " ^ Dif^') 



^ ^ {-lfcl!;px'^^'\ukhz)x^''\^khz)e'^'''-^''^'"^' 

(39) 

where /3 = ii fi = p and /3 = 1 if // = s. 

A key distinction between the "direct" and "exterior" 
methods can be illustrated through the example being 
considered here. Let us assume that the half space L is 
composed of multiple planar films. In the direct method 



35l . [36[ , the contribution of each layer is evaluated sepa- 
rately and added subsequently to determine the heat flux 
in the vacuum layer due to half space L. This requires 
finding the DGFs when r belongs to the vacuum layer and 
r lies in each of the thin films that makes up the half spaces 
L. On the other hand, the exterior method, which is used 
to derive Eq. [JHland Eq. [HI requires only knowledge of the 
DGFs when r and r belong to the vacuum layer. Though 
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Figure 2: Two multilaycrcd half spaces at Ti and T2 separated by a 
vacuum gap. 



one can derive the same results, but with more algebraie 
manipulations, for layered media, the ease of using the 
exterior method should beeome apparent when one tries 
to model near-field radiative transfer in more complicated 
geometries, for instance between two coated spheres. 

5.1.1. Two parallel half spaces: Radiative transfer 

For the two half spaces in Fig. [2l the surface normal 
vectors n(r) and n(r) in Eq. |30] are given by z and —z 
respectively. G'e(r, r), G„i(r, r), Ge{v,v), and GA/(r, r) 
derived using Eq. [33] and Eq. [31] are substituted into 
Eq. [30] to obtain, T'^'J^, the generalized transmissivity of 
energy transfer between planar half spaces (the superscript 
pp is short for "planar-planar" and is used to indicate the 
type of objects). After some manipulations, T^']^ can be 



shown to be: 



UJ I C 



27r 



/kpdkp '^-^ 
27r 



AQiRip^iR'^')e 



H=s,p 



(40) 



This expression for T^^2(j^) is in agreement with ex- 
pressions for transmissivity of energy transfer across half 
spaces . For the case of two homogeneous half spaces, we 
have also confirmed that volume integral expression (Eq. 
I16cp as well as surface integral expressions (Eq. [23 Eq. 
130)) yield the same result. Two interesting features of Eq. 
[30] need to be emphasized: (1) Eq. [3D] is valid for energy 
transfer not just between two homogeneous half spaces but 
also between two half spaces comprising planar thin films, 
(2) It is valid for isotropic materials with electric as well as 
magnetic polarizabilities, i.e., with frequency dependent e 
and ^. 

5.1.2. Two parallel half spaces: Non-equilibrium pressure 
Since we have not been able to derive a generalized 
transmissivity for momentum transfer, we use Eq. [5] and 



Eq. [9] to derive the van der Waals pressure in the vacuum 
gap. Because the film is perpendicular to the z direction, 
the van der Waals pressure is given by the zz component 
of the stress tensor, azz- Using the expressions for G'e(r, r) 
and G£;(r,r) in Eq. [TO] and Gm(r,r) and GM(r,f) in Eq. 
1201 a generalized transmissivity for momentum flux from 
L to i? is given by: 



u}/c 



kpdkp kii2 



E 

fj.=.i,p 



{l-\R\t^m + \R^^>n 



(M)|2^ 



-I- 



kpdkp |A;/j.2 
27r UJ 



■E 

tJ.=s,p 



4S(4';03f?(4^,0e-2|^'>^l' 
(41) 



The superscript m in T^J^2{^) stands for "momentum." 
Using the same notation as Antezza et. al. [s^, the 
non-equilibrium pressure in the vacuum layer due to tem- 
perature Tl of half space L while Tr = OK, denoted 
by PneqiThiTji = 0,/), is calculated using the formula 



-„e,(ri,TH = 0,0 - /|^e(c.,Ti)T;^^/(c.). 



It is in- 



teresting to note that Eq. [33 which is valid for half 
spaces with arbitrary e and coincide with the expres- 
sions for non-equilibrium van der Waals pressure and ra- 
diative transfer derived in Ref. [s^ and Ref. Q even 
though the authors of Ref. [s^ and Ref. Q derived it 
only for the case when = 1 everywhere. 

5.2. Agreement with theory of blackbody radiative transfer 
Computing radiative transfer between two arbitrarily 
shaped isotropic objects using Eq. [53 Eq. [57] or Eq. 1301 is 
computationally involved because of the need to compute 
the appropriate DGFs. However, for one class of objects, 
namely blackbodies (or objects those can be approximated 
as blackbodies), the expression for generalized transmis- 
sivity derived here can be used to obtain useful results 
irrespective of the shape. Blackbody radiative transfer is 
derived from Planck's theory of blackbody radiation and 
Kirclioff's laws, both of which are consequences of ther- 
modynamics applied to relatively simple electrodynamical 
systems (for example, photon gas in a piston with perfectly 
reflective walls). Using thermodynamic arguments, this 
idea is generalized to arbitrarily shaped objects to yield 
Qbb^ = AiFi^2crsB {T^ -T^), where asB is the Stefan- 
Boltzmann constant, the superscript bb stands for "black- 
body", and F1.2 is the view factor between objects 1 and 
2. The view factor F1.2 between the two objects in Fig. [T] 
is given by: 



Fi. 



—— J) dr J) dr 
All J 

Si S2 



(42) 



where R = (r — f)/|r — f |, R = i?R = r — r, r G 5*1, and 
r € S2, and Ai is the area of Si . 
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A blackbody is one that absorbs all radiation incident 
on it and scatters none. For an object in vacuum (eh ~ 
= 1), this can be achieved by a region of space (the 
blackbody) with permittivity and permeability given by 
£ = 1+iS, /i = 1+17 such that (5, 7 0, ensuring that there 
is no scattering by the object. The nominal dimension L 
should be such that 5kL ^ 1 or jkL ^ 1, where k = ui/c 
and c is the speed of light in vacuum, ensuring that all 
the radiation entering the object is absorbed. Because the 
properties of the objects differ infinitesimally from that of 
the host medium, scattering can effectively be neglected 
and the DGFs, Ge(r, ? ) and G,„(r, r), are simply given by 
the DGF in free space, which is: 



transmissivity between two blackbodics in the far-field is: 



Go(r,r) 



(/ - RR) 1 



khR klR^ 



)+ 



khR klR^ 



(43) 



When the spacing between objects is large compared to 
the wavelength, Eq. |43] reduces to: 



Ge(r,r) = G„,(r,r) = Go(r,r) 



AttR 



(/ - RR) (44) 



Similarly, G£;(r,r) and Gj\/(r,r) can be written as: 



V X Go(r,f) = Go(r,r) = ikh 



AttR 



(R X /) (45) 



To derive the generalized transmissivity between two 
blackbodies, the following derivations are useful: 

SRTr [(ni(r) x fih^e{r,r)) ■ (M^) x e,,^„,(f , r))*] 
=5Rrr[(ni(r) x^o(r,r)).(n2(r) x^„(r,r))*] 
[ni(r) X (T - RR)]„- [n2(r) x (7 - RR)],, 



[fii(r) 



(47ri?)2 
y[n2(r) X {§§ + 



(46) 



(47ri?)2 

_ -2(ni(r).R)(n2(r)-R) 

(47ri?)2 

5RTr [(fii (r) x^M(r,r)) • (n2(r) x ^M(r, r))*] 
= KTr [(fii(r) X %{r,r)) ■ (M^) x %{r,r)y 



[fii(r) X (R X /)]y [n2(f) x (R x /)],, 



[ni(r) X 



(47ri?)2 
J^j[n2(f) 



(47) 



(47ri?)2 
-2(ni(r).R)(n2(r).R) 



(47ri?)2 

Substituting Eq. gHand Eq. |17]in Eq. [501 the generalized 



dr J) dr 

Si S2 
L^2 



-ni(r) -Rj (^n2(r) -R 
27r2i?2 



(48) 



27rc2 



1,2 



and Q12 is given by: 



--AiFi,2'JsB{Tt-Ti) 



(49) 



The agreement of obtained results in the examples of Sec. 
15.11 with earlier published works [s^ and the above 
derivation for the radiative heat transfer between black- 
bodies attest to the correctness of our definitions of the 
generalized transmissivity (Eq. Eq. [571 and Eq. 150)) . 
Extension of the proof given here to obtain Eq. [49] to the 
case of gray body radiative transfer is not a simple one 
because we consider only specular refiection at surfaces. 

6. Summary 

In this paper, we have developed a dyadic Green's func- 
tion formalism to determine radiative heat transfer and 
non-equilibrium van der Waals/Casimir forces between ob- 
jects of arbitrary shapes, sizes, and with frequency de- 
pendent dielectric permittivity and magnetic permeabil- 
ity. The cross-spectral densities of electromagnetic fields 
are necessary to evaluate Poynting vector and electromag- 
netic stress tensor from which radiative transfer and forces 
between objects can be evaluated. Using Rytov's fiuctua- 
tional electrodynamics, expressions for cross-spectral den- 
sities in terms of volume integrals of products of DGFs are 
obtained. Green's identities for dyadic functions are then 
used to convert these volume integral expressions into sur- 
face integrals of products of tangential components of the 
DGFs on the surfaces of scatterers. The spectral radiative 
transfer between two objects is described in terms of a sin- 
gle quantity - the generalized transmissivity - which can 
be represented as a double integral of products of tangen- 
tial components over the surfaces of the two objects. The 
spectral integral of Tf_^2('^) weighted by the temperature 
derivative of the Bose-Einstein function yields the thermal 
radiative conductance between the two objects. In the 
geometric optics limit, the thermal radiative conductance 
between two blackbodies, as derived from the generalized 
transmissivity, is shown to agree with the predictions of 
the classical theory of radiative transfer. 

While many computational methods, such as finite ele- 
ment method, vector eigenfunction expansion, or T-matrix 
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method, can be employed to compute rf_^2('^)j the sur- 
face integral expression in Eq. [30] points towards the 
surface integral equation method (SIEM) as best-suited 
for the purpose. The main advantage of a surface inte- 
gral equation based method is the potential reduction in 
computational cost due to restriction of the discretization 
domain to a surface rather than a volume. This advan- 
tage becomes more important as the size of the object 
becomes larger. After the submission of this work, we 
were made aware of a recently released free software ti- 
tled SCUFF-EM (Surface CUrrcnts Field Formulation of 
Electro-Magnetism) based on a surface integral formu- 



lation of the electromagnetic scattering problem [4] 
The similarities between photon and phonon transport 
lead us to wonder whether such a surface integral formula- 
tion for phonon energy transport in mesoscale structures 
in the harmonic limit is also possible. 

This work has been funded partially by National Science 
Foundation through Grant CBET-0853723. 
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